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1 Introduction 



The exploration of the nature of the phase space that we live in, is a naturally 
attractive endeavour. On astronomical scales, this translates to an exercise 
in understanding the state of the phase space in the neighbourhood of the 
Sun. Of course, for a long time, this was impeded by the dearth of sufficiently 
large samples of proper motion data of the nearby stars. However, with the 
(— > . information available in Woolley's catalogue [I], Agris Kalnajs concluded the 

local velocity space to be basically bimodal [2] and attributed this behaviour 
to the rough proximity of the Sun to the Outer Lindblad Resonance of the 
central bar in the Galaxy (OLRb). The observational domain expanded with 
the availability of the transverse velocity information of stars in the vicinity of 
the Sun, as measured by Hipparcos. Since then many workers have attempted 
. to chart the local phase space distribution DF [31 [H [5] and investigate the 

' origin of the structure in the same [U [H] [3 [9] . 

Here we pursue the hypothesis that the observed state of the local phase 
space owes to the dynamical effect of Galactic features such as the central bar 
or the spiral pattern. However, dynamical modelling of the Solar neighbour- 
hood is incomplete without the inclusion of the effects of both the bar and 
the spiral pattern; after all, a major resonance of the bar occurs in the vicin- 
ity of the Sun, (as suggested by the Kalnajs Mechanism) and the inter-arm 
, separation in the Galactic spiral pattern falls short of the average epicyclic 

' excursion of a star at the Solar radius [TU] . 

However, this joint perturbation has been considered only occasionally 
jTUJ [8] . That too, it is contentious if the Galactic spiral arms should attach 
themselves to the ends of the bar or rotate with a pattern speed (O s ) that 
is distinct from that of the bar (i?f>). The former scenario was invoked to 
produce self-consistent models of the galaxies NGC 3992, NGC 1073, and 
NGC 1398 llj while the best-fit model for NGC 3359 was reported to manifest 
distinct pattern speeds for the bar and the spiral [12] . Such dissimilar pattern 
speeds were considered for the first time by [TU] in her modelling of the Solar 
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neighbourhood. The effect of the relative separation between the locations of 
the resonances due to the Galactic bar and the spiral pattern is expected to 
bear interesting dynamical consequences. 

In this article, we report the results of modelling of the local phase space, 
undertaken with the aim of disentangling the mystery about the origin of the 
observed phase space structure. In particular, we address the question of the 
imperativeness of the presence of chaos, in order to explain this structure. 

2 Local Phase Space 

The exploration of the physics responsible for the observed DF, follows the 
exercise of density estimation in phase space, given the observed velocity data. 
Assuming the Galactic disk to be ideally flat, we approximate the phase space 
as given by four coordinates: the two plane of the disk spatial coordinates (ra- 
dius R and azimuth <j)) and the (heliocentric) radial and transverse velocities 
(U and V, respectively). U is measured positive in the direction of the Galac- 
tic centre while V is positive along the sense of Galactic rotation. Thus, our 
simulated velocity space is 2D in nature. 

The stellar U — V distribution is noted to be highly non-linear and multi- 
modal. [TU] uses the same distribution that was used in 0]; this construction 
uses V of single stars with distance d < 100 pc and dispersion of parallax 
(j(tt) < O.Itt from the Hipparcos Catalogue and U of the 3481 stars in the 
Hipparcos Input Catalogue. This measured kinematic data, when smoothed 
by the bisymmetric adaptive kernel estimator, (as used by [5]) suggests five 
major clumps in the local velocity U — V plane (see Figure 2 in [ffi]); these 
have been identified with the Hercules, Hyades, Pleiades, Coma Berenecius 
and Sirius stellar streams or moving groups. It is the Hercules stream that 
corresponds to the smaller mode while the bigger mode constitutes the other 
four groups. 

It is accepted that other choices of the kernel might have given rise to a 
different picture of the U—V plane. This is expected to be more the case at the 
outer parts of the distribution, where there are relatively fewer stars than near 
the central parts (between -50 km s _1 < V < 10 km s _1 and -50 km s _1 < U < 
50 km s _1 ). To check if the suggestion of the five streams is an idiosyncracy 
of the density estimation procedure, other kernel estimators were tried; it was 
concluded that this did not affect the above mentioned central part of the 
distribution in question, though the outer parts were indeed affected. 

It merits mention that in the course of the density estimation procedure, as 
suggested by [5], the calculation of the local bandwidth Aj, at the i th velocity 
data point, involves achieving convergence, by iteratively toggling between 
the distribution estimate at this point and \. The effect of varying the initial 
seed for the density estimate, on the result of this iteration, has however not 
been thoroughly investigated or reported; subsequent workers [HQS] have used 
values provided by [5]. 
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Thus, accepting that the observed structure in the local velocity space 
has been robustly estimated, test particle simulations (2-D by nature) were 
carried out to constrain the model that best fit the observations. 



3 Simulations 

The simulations are test-particle in nature. This form of constricted N-body 
sampling of the stars is sufficient in the outer parts of disk galaxies and is useful 
when there is a multi-dimensional parameter space that need to be scanned. 
In the simulations, a bunch of initial phase space coordinates (^3500, i.e. 
the number of stars used to construct the local U — V diagram) are evolved 
in the potential of the Galactic disk, as perturbed by the potentials of the 
bar and/or the spiral pattern. The disk is assigned an initial phase space 
distribution function given by the doubly cut-out forms used in [13]; this 
ensures an exponential surface mass density profile. We choose to ascribe a 
Mestel potential to the disk, in order to recover a flat rotation curve; the 
warmth parameter in this potential is maintained sufficiently high to ensure 
the values of radial and transverse velocity dispersions and vertex deviations 
of stars in the vicinity of the Sun. The choice of the perturbation potential 
is motivated by their geometries: the bar is modelled as a rigidly rotating 
quadrupole while the spiral potential is approximated by a logarithmic spiral. 

3.1 Spiral Characteristics 

The spiral pattern that we use in our simulations is considered detached from 
the ends of the bar. In line with the gas dynamical model of the galaxy, [14], we 
choose the bar pattern speed to be more than double that of the spiral. This 
is corroborated by [15] who also suggests an upper-limit of 25kms~ 1 kpc _1 on 
f2 s . We experiment with three distinct ratios of S7 S : Sib - 18/55, 21/55, 25/55. 

We choose this spiral to be a 4-armed, tightly wound structure (pitch angle 
= 15°), along the lines of [16] . This implies that the ILR of this 4-armed spiral 
(ILR S ) occurs outside, on top of and inside OLRb, for the choices of Sl s : fib 
- 18/55, 21/55 and 25/55, respectively. Naturally, the 21/55 model raises 
interest, given that in this case, resonance overlap occurs, ensuring global 
chaos. [S] had invoked this scenario to explain the observed stellar streams in 
the vicinity of the Sun. 

3.2 Some Technicalities 

In our scale-free disk, all lengths are expressed in units of the co-rotation 
radius of the bar [Rcr). All pattern speeds are also expressed in terms of the 
bar pattern speed, which is set to unity. 

We choose to record orbits in an annular region between Rj Rcr — 1.7 to 
R/ Rcr = 2.3. The significance of the value of 1.7 is that the Outer Lindblad 
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Resonance of an to = 2 perturbation occurs at R/R C r = 1-7, in a background 
Mestel potential. 

All orbits are recorded in the frame that is rotating with the bar. Azimuth 
<j> = coincides with the bar major axis. In the bar+spiral simulations, where 
the figure of the spiral is not static in the frame of the bar, we record the 
orbits only when, at R — Rcr, the potential of the spiral is maximum. At the 
end of the simulations, each orbit is sampled (1000 times) in time, which is 
equivalent to sampling in phase, assuming ergodicity to be valid, to provide 
us the output. 

4 Results 

In our simulations, we work with 5 models. These are: 

• bar only model with no perturbation from the spiral: model bar — only. 

• bar and spiral, with f2 s : J? b =18/55: model bar + sp 18 . 

• bar and spiral, with S7 S : J?{,=21/55: model bar + sp^i- 

• bar and spiral, with fi s : J?b=25/55: model bars + P25- 

• spiral only model with no perturbation from the bar: model spiral — only 

The output orbits are first put on a regular polar grid; each of our R — <f> 
bins are chosen to be 0.025Rcr wide in radius and 10° wide in azimuth. At 
each relevant R — <j> cell, the recorded velocities are put on a regular Carte- 
sian U — V grid. With this stored kinematic data, velocity distributions are 
prepared at each R — <f> bin. When the simulated velocity distribution at a 
given disk location {f s {R, </>))is found to match the observed distribution (f Q ) 
well, we refer to such a velocity distribution as a "good" distribution and the 
corresponding R — <j> bin is called a "good" location. 

The quality of overlap between f Q and f s (R, (f>) is quantified by a goodness- 
of-fit technique (discussed below). All this is undertaken individually for each 
of the 5 models. 

4.1 Goodness of Fit Technique 

We test the null hypothesis: the observed U — V data is drawn from f s (R, (f>) . 
This testing is done by estimating the p-value of a test statistic S. p-value 
measures the probability of how unlikely it is that a null hypothesis is true 
by fluke; thus, low p-values indicate bad fits to the data. Though extensive 
literature exists to suggest various options for S, we find that the definition of 

5 as the reciprocal of the likelihood L, suffices in our case. Here, the likelihood 
of data D, given the simulated distribution f s (R, (f>) is: 

L{D\f s )=]\v R ,^U u V 3 ), (1) 

where VR^(Ui, Vj) is the value of f s (R, </>), in the (i,j) th U — V bin and the 
product is performed over all the U and V velocity bins. We decide to choose 
the pre-set significance level of 5%, for the acceptance of the null hypothesis. 
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4.2 Identification of "Good" Locations 

For a given model, the quality of overlap between f a and f s (i?, <fi) is given 
by the p-value of the test statistic in that R — <fi bin. This is referred to 
as p — value(R, cf>). Distributions of p — value(R, <p) are shown in Figure 1, 
for models bar — only and spiral — only. It is to be noted that the function 
p~value(R, <j)) for the resonance overlap (bar + sp2i) and spiral — only models 
is distinct from that for the other three models. This is attributed to the chaos 
inducing efficacy of the spiral perturbation. 




From the distribution p — value(R, </>) , "good" locations are identified as 
those R — <f> addresses, at which the p-value attains the maximum, i.e. 100%. 
This distribution of the "good" locations can be marginalised over <j> to allow 
us to quantify the ±1-<t range that defines the best radial locations. Similarly, 
marginalising over R gives the constraints on the best azimuthal location for 
the observer, in a frame where </>=0 implies bar major axis. Of course, these 
are then constraints on our location, i.e. the Solar position, in the frame of the 
bar. Such marginalised distributions (with R and </>) are shown in Figure 2, 
for model bar + sp25 ■ This tells us that 

• the best radial location for the observer is R/Rcn— 2.0875(j'25- Equating 
the Solar radius (of 8 kpc) with the medial value, we scale the bar pattern 
speed to 57. 4^1 km s _1 kpc -1 (using circular velocity at Sun=220 km s _1 ). 

• the best fit values for the bar angle is [0°, 43°], with a medial value of 6°. 

We find that all the models that include the influence of the bar and the 
spiral pattern are able to produce the observed stellar streams in the vicinity 
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Radius (R CR ) Azimuth (degrees) 

Fig. 2. The distribution of the maximal p-value locations is marginalised, over 
azimuth and radius, to yield best-fit values of observer radial and azimuthal location 
(respectively), in order for the observer's local velocity distribution to best match 
the observed distribution. 

of the Sun, at different disk locations. Though the bar — only model is also 
capable of the same, it has to be discarded in light of the dynamical reasoning 
that no model of the Solar neighbourhood is complete without taking the effect 
of the spiral into account ( I10j). The spiral — only model is rejected since it 
fails to warm up the local patch in the disk enough, to guarantee velocity 
dispersions of the order of those observed in the Solar neighbourhood. 

5 Origin of the Splitting of the Larger Mode 

On the basis of the results above, we conclude that the splitting of the larger 
mode in the local U — V distribution is not due to resonance overlap since 
bar -(-spiral models that do not impose resonance overlap are also successful 
in producing the observed phase space structure. On the other hand, this 
splitting is suggested to be due to the interaction of chaotic orbits and orbits 
that belong to families corresponding to the minor resonances due to the spiral 
and the bar (such as the outer 1 : 1 resonance of the bar) . 

As for the overall bimodal character of the local phase space, a rudimentary 
orbital analysis confirms that in line with the Kalnajs Mechanism, the overall 
bimodality is indeed due to scattering off OLRf,; the Hercules stream appears 
to be built of quasi-periodic orbits belonging to the anti-aligned family while 
the anti-aligned family was not spotted in the other larger mode. (The other 
mode manifests orbits from the aligned family). 
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Fig. 3. Surface of section of orbits at effective energy value of -1, from the bar + spig 
model. The strongly chaotic orbits are in medium grey shading, while the regular 
orbits are in the darkest shade. The intermediate shading corresponds to the weakly 
chaotic orbits. 



6 The Role of Chaos 



[4] has suggested bar-induced chaos as the cause for the local streams. How- 
ever, in a chaos quantification exercise that is currently underway (Chakrabarty 
& Sideris, submitted to A&cA), the orbits in the bar — only model are found 
to be regular; it is noted that it is the presence of the spiral that triggers 
the onset of chaos. Figure 3 indicates the relative fraction of chaotic orbits, 
over regular and weakly chaotic orbits, in the bar + spis models. Similar sur- 
faces of section for the bar — only models, at even higher energies indicate no 
irregularity. 

It is of course possible that a stronger bar would induce greater irregularity 
than in our work, but the point is that the local phase space structure can be 
emulated with the bar strength used herein. 



7 Summary & Future Work 



In this contribution, we have attempted to understand the state of the local 
phase space via dynamical modelling that includes the influence of the bar 
and the spiral pattern in the Galaxy. The clumpy structure of the 2-D local 
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velocity distribution was attributed to the major and minor resonances of 
these perturbations, as well as the emergence of chaoticity triggered primarily 
by the spiral potential. Only the models that take the gravitational potentials 
of both bar and spiral are found to offer viable representations of the Solar 
neighbourhood dynamics; overlap of the results from the successful models 
indicate a fast bar, with a pattern speed of 57.4+ 2 ' 5 km s _1 kpc _1 and a bar 
angle that lies in the range [0°, 30°]. 

The currently used p- value formalism is not satisfactory on account of the 
lack of automativeness and inability to rank or grade the models. In fact, we 
are hoping to improve upon this technique, with a Bayesian measure of the 
quality of a simulation. This is to aid the expansion of the relevant parameter 
space, particularly, in the inclusion of the effects of the halo and progression 
to a full six dimensional phase space. 
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